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1 Abstract 


The EPOXI Discovery Mission of Opportunity reused the Deep Impact flyby spacecraft 
to obtain spatially and temporally resolved visible photometric and moderate resolution 
near-infrared (NIR) spectroscopic observations of Earth. These remote observations pro- 
vide a rigorous validation of whole disk Earth model simulations used to better under- 
stand remotely detectable extrasolar planet characteristics. We have used these data to 
upgrade, correct, and validate the NASA Astrobiology Institute’s Virtual Planetary Labora- 
tory three-dimensional line-by-line, multiple-scattering spectral Earth model (Tinetti et al., 
2006a, b). This comprehensive model now includes specular reflectance from the ocean and 
explicitly includes atmospheric effects such as Rayleigh scattering, gas absorption, and tem- 
perature structure. We have used this model to generate spatially and temporally resolved 
synthetic spectra and images of Earth for the dates of EPOXI observation. Model param- 
eters were varied to yield an optimum fit to the data. We found that a minimum spatial 
resolution of ~100 pixels on the visible disk, and four categories of water clouds, which 
were defined using observed cloud positions and optical thicknesses, were needed to yield 
acceptable fits. The validated model provides a simultaneous fit to the Earth’s lightcurve, 
absolute brightness, and spectral data, with a root-mean-square error of typically less than 
3% for the multiwavelength lightcurves, and residuals of ~10% for the absolute brightness 
throughout the visible and NIR spectral range. We extend our validation into the mid- 
infrared by comparing the model to high spectral resolution observations of Earth from 
the Atmospheric Infrared Sounder, obtaining a fit with residuals of ~7%, and brightness 
temperature errors of less than IK in the atmospheric window. For the purpose of under- 
standing the observable characteristics of the distant Earth at arbitrary viewing geometry 
and observing cadence, our validated forward model can be used to simulate Earth’s time 
dependent brightness and spectral properties for wavelengths from the far ultraviolet to the 
far infrared. 
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2 Introduction 


After an initial decade dominated by the discovery of Jupiter-mass planets, the next 
frontier of exoplanet research will be the detection and characterization of terrestrial- 
mass planets. Within the next few years, NASA’s Kepler spacecraft will make the first 
comprehensive estimates of the prevalence and nature of extrasolar terrestrial planets 
(Borucki et al., 2003), while searching for Earth- mass planets in the habitable zones 
of their parent stars (Basri et ah, 2005). In the coming decades more ambitious 
planet detection and characterization missions for habitable Earth-mass planets are 
planned, such as NASA’s Terrestrial Planet Finder mission (Beichman et ah, 1999). 
These missions will be designed to detect and characterize nearby habitable planets, 
with the capability to obtain direct imaging, and photometric and spectroscopic data 
for extrasolar terrestrial planets 

The observational challenges inherent in characterizing a terrestrial exoplanet are 
significant, and carefully considered trade offs must be made to maximize the sci- 
ence return. Even with the most ambitious telescopes planned, terrestrial exoplanets 
will remain faint, spatially unresolved point sources. The principal challenge is to 
determine the minimum and optimum sets of observational parameters that can best 
characterize the environment of an unresolved planet, which may be spatially inho- 
mogeneous, cloud covered, and temporally variable. For example, the combination 
of temporal resolution and multi-wavelength photometry could disentangle phase or 
rotation dependent differences in surface properties from variable cloud cover. The 
resulting maps could discriminate between large scale surface inhomogeneities such 
as continents and oceans (Palle et ah, 2008; Cowan et ah, 2009; Kawahara and Fu- 
jii, 2010; Fujii et ah, 2011). Disk integrated spectroscopy can potentially determine 
globally averaged atmospheric and surface composition to verify habitability and to 
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search for global evidence of life in the planetary environment (Seager et ah, 2005; 
Montanes- Rodriguez et al., 2006; Meadows, 2006). 

New tools are needed to obtain quantitative information about the environments 
of terrestrial planets that can only be studied as unresolved point sources. A typical 
approach to understanding a world from disk-integrated observations consists of a 
“forward model”, an “instrument model”, and an “inverse model.” The forward model 
is typically a radiative transfer model designed to generate a synthetic spectrum, given 
an assumed surface-atmospheric state and viewing geometry. The instrument model 
simulates the spectral and spatial resolution and other properties of the observing 
system. The inverse model adjusts surface or atmospheric state to yield a better 
fit to the observations. Given a candidate observing system design, refinements in 
both forward models and inverse methods are needed to fully exploit the information 
content of disk- integrated observations of terrestrial planets. Most existing forward 
models are designed to analyze soundings taken with adequate spatial resolution to 
yield spatially homogeneous sounding footprints. Forward models designed for surface 
or “clear sky” remote sensing applications rarely perform well with cloudy soundings. 
Those designed for land remote sensing observations rarely simulate the reflection 
from the ocean surface. In short, few if any forward models have been designed to 
yield accurate observations over the full range of solar illumination angles, observation 
angles, or surface and atmospheric properties present in a single, integrated, full disk 
observation of an extrasolar planet. Here, we describe a forward model designed 
specifically to simulate disk- integrated observations of Earth, which can serve as a 
theoretical “laboratory” for the accurate simulation of Earth’s appearance at arbitrary 
viewing geometries and phases. These simulations can be used to explore and identify 
the best conditions under which to search for planetary characteristics of habitability 
and life, such as the presence of an ocean or a photosynthetic biosphere, and can 
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also be used to generate test data to challenge proposed observational and retrieval 
techniques for extrasolar planet characterization. 

Remote sensing observations provide a glimpse of Earth’s appearance from space. 
Our planet is well studied by an armada of satellites, covering several wavelength re- 
gions from the ultraviolet through microwave with high temporal and spatial sampling 
of Earth’s photometry and spectroscopy (Hearty et al., 2009). However, these data 
sets are often not ideal for studying Earth as an astronomical target because of the 
sheer data volume that must be manipulated to produce a global, or disk-integrated, 
view, and because most satellites are in Sun-synchronous low Earth orbits, which 
view the Earth over a limited range of viewing geometries and times of day. By com- 
parison, spatially and/or spectrally resolved forward models can readily simulate the 
full disk Earth from observing geometries, solar illuminations, and wavelengths that 
may not be accessible to existing Earth observing satellites. 

Existing Earth models for exoplanet characterization studies are largely dominated 
by computationally inexpensive specular reflectance models (e.g., Ford et ah, 2001; 
Williams and Gaidos, 2008). These models do not include atmospheric absorption 
and scattering, and are most effective at visible light wavelengths, where they can be 
used to model the photometric variability of the rotating Earth, and to determine the 
detectability as a function of phase of specular reflection or “glint” from the ocean sur- 
face. Similar models (Pallc et ah, 2003; Oakley and Cash, 2009) utilize bi-directional 
reflectance functions that are designed to match data which have been measured by 
Earth observing satellites (e.g., Manalo-Smith et ah, 1998). Spectral Earth mod- 
els that include absorption by some atmospheric species and which simulate Earth 
spectra using weighted averages of independent, one- dimensional component spectra 
representing cloudy and clear sky scenes have also been developed (Woolf et ah, 2002; 
Turnbull et ah, 2006; Stam, 2008). Montanes- Rodriguez et ah (2006) modeled Earth’s 
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spectrum using a line-by-line radiative transfer model which was limited to simulating 
Earth’s environment in one dimension, using globally averaged atmospheric, surface, 
and cloud properties. In their model, standard atmospheric composition and tempera- 
ture profiles were assumed and Earth’s spectrum was obtained by averaging different 
component spectra based on data from the International Satellite Cloud Climatol- 
ogy Project. The primary limitation associated with one- dimensional approaches to 
modeling Earth’s disk integrated spectrum is that these models do not capture lat- 
itudinal and longitudinal variations in the composition and temperature of Earth’s 
atmosphere, and thus cannot be used to quantify the impact of spatial variations in 
temperature and composition on the information content of simulated observations. 

The EPOXI mission, a reuse of the Deep Impact flyby spacecraft, has recently 
provided a very rigorous data set for validating full disk forward models of Earth’s 
spectrum prior to predictive use. The EPOXI observations of the distant Earth 
(Livengood et ah, 2011) provide several days of time resolved, multi- wavelength vis- 
ible photometry and near-infrared (NIR) spectroscopy of Earth’s disk through a full 
24-hour rotation period. Previously, disk integrated spectra and photometry of Earth 
had been obtained from limited, often single, measurements of Earthshine and from 
planetary spacecraft. Earthshine is reflected light from Earth illuminating the night 
side of the Moon (Goode et ah, 2001; Arnold et ah, 2002; Woolf et ah, 2002; Turnbull 
et ah, 2006; Montahes- Rodriguez et ah, 2006), and it is constrained by viewing geom- 
etry to restricted phases and temporal durations. Snapshot or flyby observations of 
the Earth have also been obtained from spacecraft en route to other planets, for ex- 
ample from the Galileo spacecraft (Sagan et ah, 1993) or from Mars Global Surveyor 
(Christensen and Pearl, 1997). 

Using the new EPOXI data set, we validate and present the most comprehensive 
spectral Earth model to date for the prediction of the photometric and spectroscopic 
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characteristics of Earthlike exoplanets. This model is a forward model, which is used 
to simulate the appearance of Earth to an observer, for the purpose of exploring the 
detectability of Earth’s planetary characteristics as a function of observational geom- 
etry and time. Forward models, such as the model presented in this work, are distinct 
from, but complementary to, retrieval models designed to retrieve atmospheric char- 
acteristics from observations of extrasolar planets (e.g., Madhusudhan and Seager, 
2009). In particular, it is important that forward models be as realistic as possible so 
that they accurately represent the appearance of planet, and, as a result, do not bias 
the observed planetary properties inferred when using the model as a predictive tool. 

A previous, more limited version of this model, published in Tinetti et al. (2006a, b), 
has been corrected, and significantly updated and improved to allow accurate predic- 
tions of Earth’s time dependent photometric and spectroscopic brightness, on hourly 
to yearly timescales, through realistic modeling of the radiative effects of a surface 
ocean, atmosphere, and clouds. The previous model allowed for an arbitrary scaling 
of its input cloud coverage data, and it also used the optical thickness of clouds as 
free parameters. By tuning the previous model, snapshot observations of Earth could 
be reproduced and the model could then be used to explore certain characteristics of 
Earth, such as how the planet’s brightness changes with phase. Time resolved data 
were not used for validation. 

In the process of validating the model against the time and phase resolved EPOXI 
observations, we have significantly upgraded the model to self-consistently utilize 
satellite derived cloud data. Cloud coverage is now taken from observations by Earth 
observing satellites and is no longer scaled in order to reproduce observations. The 
optical thickness of clouds in our new model is also provided by Earth observing 
satellites, rather than being tunable free parameters. We have also corrected an error 
in the model presented in Tinetti et al. (2006a, b) which effectively confused forward 
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scattering with backward scattering, thus causing the model to produce unphysical 
simulated observations of Earth. Our Earth model is based on a fully multiple- 
scattering, line-by-line radiative transfer model, SMART (Meadows and Crisp, 1996; 
Crisp, 1997), which is at the core of the exoplanet simulations generated by the 
NASA Astrobiology Institute’s Virtual Planetary Laboratory (VPL). The new, self- 
consistent treatment of clouds in our model has allowed us to match the EPOXI 
Earth observations, which span almost three months in time and a variety of phases, 
without tuning from one dataset to the next (for validation against Earthshine data 
over a wider range of phases, see Robinson et ah (2010)). Our validated model is 
capable of simulating Earth’s wavelength dependent temporal variability and absolute 
brightness, for any given viewing geometry and phase, over wavelengths from the 
ultraviolet to the far infrared (IR). 

3 Description of EPOXI Earth Observations 

The EPOXI Earth data sets used for our model validation were acquired with the 
High Resolution Instrument (HRI) on board the Deep Impact flyby spacecraft, and 
are described in (Livengood et ah, 2011). The HRI is a 0.3m // 35 telescope with nine 
square-bandpass filters and a NIR spectrometer (Hampton et ah, 2005). Seven filters 
are ~ 100 nm wide and are centered at 350, 450, 550, 650, 750, 850 and 950 rim, 
segmenting the visual spectral range. The NIR spectrometer covers the wavelength 
range from 1.05-4.5 /im with a maximum resolution (R = A/AA) of 750 at the shortest 
wavelengths, a minimum resolution of 200 at about 2.6 /mi and a resolution of 350 
at the longest wavelengths (Klaasen et ah, 2008). 

Earth was observed with the HRI in three separate 24-hour periods: 2008-Mar- 18 
18:18 UT - 2008-Mar-19 18:18 UT, 2008-May-28 20:05 UT - 2008-May-29 20:05 UT, 
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and 2008-Jun-4 16:57 UT - 2008-Jun-5 16:57 UT. The star-planet-telescope angle 
(phase angle) for the three sets of observations were 57.7°, 75.1°, and 76.6°, respec- 
tively. NIR spectra of Earth were taken every two hours (30° rotation) a total of 
13 times, to cover a full rotation of Earth from the viewpoint of the spacecraft. The 
NIR spectrometer slit was oriented perpendicular to the terminator, due to spacecraft 
constraints. To avoid saturation and to obtain globally averaged spectra, rapid scans 
were made alternating between south to north and north to south across Earth to 
limit the amount of time any part of Earth was within the slit of the spectrometer. 
Photometry in the 450, 550, 650 and 850 nm filters was taken every fifteen minutes 
(~ 4° rotation) and every hour (15° rotation) for the 350, 750 and 950 nm filters. 
Note that a subset of the March and June data sets were published by Cowan et al. 
(2009), who treated the data in a retrieval sense and performed a principal component 
analysis of the lightcurves to map the latitudinally averaged distribution of land and 
ocean on Earth. 


4 Description of Model 

The comprehensive Earth model presented in this paper uses input data from Earth 
observing satellites and a radiative transfer model to create spatially- and temporally- 
resolved spectra of Earth. Our model is a new and improved version of the VPL’s 
spectral Earth model (Tinetti et al., 2006a, b). The model uses date specific informa- 
tion on spatially resolved atmospheric composition and cloud and surface properties 
to derive a spatially and spectrally resolved datacube. This datacube can be used 
to generate both images of Earth over a given wavelength range as well as spatially 
resolved or disk integrated spectra. 

In our new model, the surface of Earth is divided into a number of pixels of equal 
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area according to the Hierarchical Equal Area isoLatitude Pixelization (HEALPix) 
model (http://healpix.jpl.nasa.gov) (Gorski et ah, 2005). HEALPix partitions a 
sphere into a number of equal area pixels, which is ideal for a planetary model where 
the exoplanet may be viewed at arbitrary viewing geometries. The resolution of the 
HEALPix model is defined by an integer N and the number of pixels used to cover a 
sphere is then equal to 121V 2 , which is a behavior inherent to the HEALPix scheme. 
A set of surface pixels are nested beneath a set of atmospheric pixels, which are 
also defined according to the HEALPix scheme. Our new model allows the surface 
resolution and the atmospheric resolution to vary independently, providing improved 
surface resolution at minimum computational cost. 

For each surface type (see Sec. 4.1) lying beneath each atmospheric pixel, we 
run a one-dimensional, line-by-line radiative transfer model, the Spectral Mapping 
Atmospheric Radiative Transfer (SMART) model developed by D. Crisp (Meadows 
and Crisp, 1996; Crisp, 1997), over a grid of solar zenith angles and observer zenith 
and azimuth angles, and for wavelengths from the ultraviolet to the far IR, to create 
a look-up table of synthetic spectra. For a given sub-solar and sub-observer latitude 
and longitude, spectra from each atmospheric pixel are taken from the look-up table, 
interpolated over solar and observer azimuth and zenith angles, and combined at the 
surface resolution of the model to create a three-dimensional spectral map of Earth. 
The spectral map can be integrated over a given wavelength range and used to create 
images of Earth, as shown in Fig. 1, or the three-dimensional spectral map can be 
integrated over solid angle to create a disk-integrated spectrum of Earth. Note that 
SMART is the same radiative transfer model used as the core to the Earth model 
presented in Tinetti et ah (2006a, b). The following subsections describe the input 
data to our Earth model. 
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4.1 Surface Properties 

We determine the spatial distribution of Earth’s surface types from a yearly averaged 
map obtained from the Moderate Resolution Imaging Spectroradiometer (MODIS) 
instruments (Salomonson et ah, 1989) aboard NASA’s Terra and Aqua satellites 
(http://www-modis.bu.edu/landcover/). In general, the largest seasonal variability 
in surface reflectivity will be due to snowfall or advancing/retreating sea ice, so the 
yearly averaged data is tailored to a specific date of observation by using eight-day- 
averaged snow cover (Hall et ah, 1995) and sea ice cover (Riggs et ah, 1999) from 
MODIS observations (http://modis-snow-ice.gsfc.nasa.gov/). 

The wavelength dependent surface reflectivity was characterized by five surface 
types: ocean, forest, grassland, desert and snow. Specular reflectance from liquid 
water surfaces in our model is simulated using the Cox-Munk glint model (Cox and 
Munk, 1954), which allows for the calculation of the bidirectional reflectance distri- 
bution function of a wave covered ocean given wind speed and direction, which are pro- 
vided by the QuikSCAT satellite (http:/ /winds.jpl.nasa.gov/missions/quikscat/index.cfm). 
Other surfaces are assumed to be Lambertian and reflect isotropically. Wavelength 
dependent albedos for non-ocean surfaces are taken from the USGS Digital Spectral 
Library (http://speclab.cr.usgs.gov/spectral-lib.html) and the ASTER Spectral Li- 
brary (http://speclib.jpl.nasa.gov/). The Cox-Munk formalism is an improvement 
in our new model over the model presented in Tinetti et ah (2006a, b), which as- 
sumed that oceans scatter as a Lambertian surfaces with an albedo near 5-6% at 
most wavelengths. 
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4.2 Atmospheric Properties 

To simulate molecular absorption in Earth’s visible and NIR disk integrated spec- 
tra and to accurately predict planetary brightness temperature in the mid-IR we 
require the three-dimensional distribution of atmospheric gases and temperatures as 
input to the model. The model includes both Rayleigh scattering by air molecules 
as well as absorption from H 2 0, C0 2 , 0 3 , N 2 0, CO, CH 4 and 0 2 . Spatially re- 
solved mixing ratio profiles for atmospheric gases are obtained from the Microwave 
Limb Sounder (MLS) (Waters et ah, 2006), the Tropospheric Emission Spectrometer 
(TES) (Beer et ah, 2001) (both aboard NASA’s Aura satellite) and from the Atmo- 
spheric Infrared Sounder (AIRS) (Aumann et ah, 2003) aboard NASA’s Aqua satellite 
(http://mls.jpl.nasa.gov/data/, http://tes.jpl.nasa.gov/data/, http:/ /airs. jpl.nasa.gov/). 

An abridged list of the species investigated by these instruments and the valid ranges 
for profile retrievals are shown in Table I (Livesey et al., 2007; Payne et al., 2009). 

Data from Aqua/AIRS and Aura/MLS are combined to produce spatially resolved 
temperature profiles. These atmospheric properties are averaged over each atmo- 
spheric pixel and resolved into 40 vertical layers prior to input to our one-dimensional 
radiative transfer model. Absorption cross sections for gases are generated using the 
HITRAN 2004 line list database (Rothman et al., 2005) (http:/ /www. cfa.harvard.edu/hitran/). 
Line profiles are simulated using a line-by-line absorption coefficient model (LBLABC) 
developed by D. Crisp (Meadows and Crisp, 1996). 

4.3 Clouds 

The reflectivities, optical depths, and the spatial and vertical distribution of clouds 
have a profound effect on Earth’s time variable spectrum. In our model, the spa- 
tial distribution of clouds is straightforwardly obtained from cloud coverage maps 
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provided by the MODIS instruments (http://modis-atmos.gsfc.nasa.gov/). MODIS 
provides only cloud phase assignments (liquid, ice and undetermined), optical depth 
measurements, and cloud top pressure. MODIS does not directly report cloud al- 
titude or reflectivity. Other data sets (e.g., CloudSat or the International Satellite 
Cloud Climatology Project) can provide more detailed information about cloud dis- 
tribution but suffer from either poor spatial coverage or a large lag time between 
data acquisition and release, which makes these data sets a poor choice for the date 
dependent simulations presented in this work. 

Not all locations on the planet have MODIS data recorded within a given 24 hour 
period. We interpolate cloud coverage and optical depths to the times of EPOXI ob- 
servations using the closest Aqua/MODIS and Terra/MODIS data that bracket the 
observation times (Wolfe, 2006). The spectral model has been improved to allow for 
an arbitrary number of cloud varieties, where a cloud variety is defined according to 
its phase (liquid or ice), altitude and optical thickness, whereas the previous model 
(Tinetti et ah, 2006a, b) used a fixed number of cloud varieties. The previous model 
used simulated as well as satellite observed cloud locations, but allowed the global 
coverage and optical thickness of each cloud type to be a free parameter. Our new 
model now Exes cloud coverage and optical depth based purely on satellite measure- 
ments, instead of allowing these to be free parameters that can be tuned to fit an 
observation. 

To model the clouds, we used the MODIS cloud phase identification to catego- 
rize a cloud at a given spatial location as either liquid or ice. MODIS has a third 
category, “undetermined” , which covers ~ 1% of the planet and is neglected in our 
simulations. Model ice clouds were placed near 8.5km (0.331 bars) altitude, and model 
water clouds were placed near 1.5km (0.847 bars). The average cloud top pressure 
in our model agrees with the average cloud top pressure reported by MODIS. Liquid 
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and ice clouds were then sub-divided using MODIS optical depth measurements into 
two sub-categories, for a total of four cloud varieties. These sub-categories allow us to 
resolve different cloud thicknesses in the model and can be thought of as moderately 
optically thick and optically thick cloud categories. The sub-division occurs at an 
optical depth of 10 for both liquid and ice clouds, since roughly half of all liquid and 
ice clouds within the MODIS data have optical depths smaller than this value. A 
weighted average of the cloud data within these sub-categories yields best-fit optical 
depths of 5 and 15 for our two liquid cloud types, and optical depths of 5 and 20 for 
our ice clouds types. This approach and the derived values held for all three dates of 
EPOXI observations. Wavelength dependent optical properties for liquid clouds were 
derived using a Mie theory model (Crisp, 1997) and ice clouds were parametrized 
using geometric optics (Muinonen et ah, 1989). The cloud scattering properties were 
assumed to be dehned strictly by cloud phase (he., liquid or ice). 

4.4 Comparison to Previous Model 

To demonstrate the significant improvements made in our new model, we compare 
spectral, phase dependent results from our model to both the model described in 
Tinetti et al. (2006a) as well as to EPOXI observations in Fig. 2. In this figure, 
we show Earth’s wavelength dependent reflectivity at visible and NIR wavelengths 
for a variety of different phases (full, gibbous, half illuminated or quadrature, and 
crescent). Both models assume realistic clouds. In general, the model described 
in Tinetti et ah (2006a) is significantly brighter than our new model, with the most 
extreme discrepancy occurring at 0.5 //m in the crescent view, where the Tinetti et. al 
model is about 400% brighter than our new model. Also, the old model is distinctly 
bluer than our new model at phases between half illumination and crescent. EPOXI 
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observations taken at gibbous phase and near half illumination (phase angles of 57.7° 
and 76.6°, respectively) and our model of the observations are shown as dashed lines 
in the gibbous and quadrature plots, demonstrating our ability to reproduce Earth’s 
brightness and reflectivity at these phases. Further validation of our model against 
measurements of Earth’s phase dependent brightness and reflectivity from Earthshine 
observations can be found in Robinson et al. (2010). 

Some of the differences between the two models shown in Fig. 2 can be accounted 
for by an error in the model described in Tinetti et al. (2006a) which rotated the angle 
between the Sun and the observer by 180°, effectively confusing forward scattering 
with backward scattering. Inserting this error into our model allows us to reproduce 
the excess brightness in the old model near full phase, and also allows us to reproduce 
the incorrect blue nature of the old model near quadrature and crescent phases. We 
have not, however, been able to reproduce the excess brightness in the old model 
at phases away from full phase, indicating that other problems exist in the model 
described in Tinetti et al. (2006a). It is important to note that the confusion between 
forward and backward scattering in the Tinetti et al. model compromises the phase- 
dependent, reflected light brightnesses and reflectivities presented in Tinetti et al. 
(2006b). 

5 Results 

5.1 Validation with Lightcurves and Spectra 

Comparisons between EPOXI observations and our model for all three epochs of 
observation are shown in Fig. 3. These plots demonstrate the brightness of Earth 
through seven EPOXI bandpasses over the duration of the 24 hour observation. Both 
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the data and the model have been normalized to their respective 24 hour averages 
through these filters. The vertical gray line marks the beginning and end of the 
observations. The discontinuity in normalized intensity between beginning and end 
is real, and is due to Earth’s time varying cloud formations, which typically produce 
a 3-5% change in the lightcurve over all wavelengths (relative errors for an individual 
bandpass are typically much less than 1%). The model generally reproduces the sign 
and magnitude of this discontinuity. The shaded region of the lightcurves for the 
May set of observations marks a Lunar transit of Earth’s disk, which is an effect not 
included in our model. The comparison between the 24 hour averaged radiance from 
the data and our model for the March and June epochs of observation is shown in 
Fig. 4. Table If shows the percent difference between the 24 hour average radiance 
data and the standard March, May, and June models (models “a”, “b”, and “c”, 
respectively). The model reproduces the data on an absolute scale of radiance to 
within the ~ 10% uncertainty in EPOXl/HRI calibration accuracy (Klaasen et ah, 
2008). Root-mean-square (RMS) errors for the lightcurve comparisons are also shown 
in Table 11. These errors measure the goodness of fit of the model lightcurves and are 
computed by comparing the data to the model for each bandpass at each observation 
within the 24 hour observing sequence. Our model reproduces the time dependent 
variability to within 3% in most cases. Note the trend in the data of larger peak- 
to-trough variability with increasing wavelength, with roughly 10-15% variability at 
shorter wavelengths and over 20% variability (in some cases) at longer wavelengths. 

In addition to fitting lightcurve time variability and absolute brightness in the 
visible, our model simultaneously reproduces the EPOXI NIR spectral data, shown 
in Fig. 5. This figure shows comparisons between our model and the EPOXI data 
for a variety of viewing geometries from the March set of observations as well as 
for data from the May and June sets of observations, which are dimmer in reflected 
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light clue to Earth phase. Residuals from the data-model comparison are also shown 
and are typically less than about 15%, demonstrating the the ability of the model 
to reproduce spectral observations on timescales from hours to months. Notable 
disagreements occur near the center of the 1.4 // m and 1.9 /jrn water bands and near 
4.1 fin l, in the short wavelength wing of the 4.3 /mi CO 2 feature. The absolute 
magnitude of the 4.1 /mi defect is much smaller than the 1.4 /jrn and 1.9 fin 1 defects. 
Instrument calibration uncertainties are typically 10%, and tend to increase below 
2.0 /m 1 and above 4.3 /jrn (Klaasen et ah, 2008). 

In Fig. 6 we show a comparison between the 24 hour average mid-IR spectrum 
from our March model and a mid-IR spectrum of Earth generated from AIRS obser- 
vations taken over the same timeframe (Hearty et ah, 2009). The AIRS instrument 
does not record full disk observations of Earth, so the AIRS spectrum is assembled 
from multiple scenes recorded over the observational period. The fact that many 
observations must be combined to produce a disk integrated spectrum implies that 
this technique cannot provide the same time resolution as the EPOXI data set. In 
general, the model reproduces the AIRS observations quite well, with residuals being, 
on average, about 7%, and with brightness temperature errors less than IK in the 
atmospheric window. 

For comparison we ran the original VPL Earth model (Tinetti et ah, 2006a, b), 
which was only validated with single snapshots of Earth, using input cloud coverage 
maps as well as atmospheric composition and temperature data for the appropriate 
dates of EPOXI observation. As the core radiative transfer model is the same in both 
Earth models, this experiment primarily tested the cloud parametrization in the two 
models. The comparisons between EPOXI data from March and the 2006 model are 
shown in Fig. 7 (model “d”). Discrepancies in the 24 hour average radiance and the 
RMS errors for the lightcnrves are also shown in Table II. The inability of the model 
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to reproduce the lightcurves indicated that a new cloud parametrization was needed, 
and demonstrated the importance of validating spectral Earth models against time 
and phase resolved data. 


5.2 Model Sensitivity 

Onr selected atmospheric and surface resolution (48 pixels and 192 pixels, respec- 
tively) and onr set of four cloud categories represent onr standard model. The follow- 
ing set of investigations aim to determine the level of detail required in the model to 
reproduce both the visible and NIR EPOXI observations while remaining consistent 
with the input data. The parameters used in these studies and a summary of the 
results are shown in Table II. Lightcnrves for a subset of the studies through three 
EPOXI filters are shown in Fig. 7. Results are presented and discussed in greater 
depth below. 

5.2.1 Cloud Categories: The Importance of Spectra 

The characteristics that define onr selected clond sub-categories are based on MODIS 
data, but the number of cloud varieties in onr model is arbitrary. While onr standard 
model utilizes four cloud sub-categories, it is useful to know how sensitive the model 
is to the chosen number of cloud sub-categories. To test this, the model was run 
with only a single clond sub-category (model “e” in Table II) and with two cloud 
sub-categories. The single cloud model uses a cloud extinction optical depth of 10 
and places the cloud in the middle of the troposphere while the two cloud model 
uses liquid water clouds with an extinction optical depth of 5 and ice clouds with 
an extinction optical depth of 15, placed at the same altitudes as the liquid and ice 
water clouds in the standard model. In all cases the characteristics of the cloud sub- 
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categories were derived from MODIS observations. Example NIR spectra from these 
low “cloud resolution” models are compared to the EPOXI data and our standard 
model in Fig. 8. In general, the models with less than four cloud sub-categories are 
poorer fits to the NIR data. The single cloud model is too bright in the continuum 
regions, underestimates water vapor absorption near 1.4 /mi, and overestimates water 
vapor absorption near 1.1 /mi. Furthermore, this model cannot reproduce the shape 
of the continuum region near 1.6 /mi, where water ice absorbs. Thus, the EPOXI 
disk integrated NIR spectra contain evidence for the presence of both liquid and ice 
water clouds. The two cloud model is an improved fit, but struggles to reproduce the 
measured intensity within the 1.1 /mi and 1.4 /mi water features. A six cloud model 
was run (not shown here) and did not offer significant improvements over the four 
cloud model when compared to the EPOXI data set. 

The visible EPOXI lightcurves are less of a challenge to models due to their 
broadband nature and the fact that these data are relatively insensitive to atmospheric 
absorbers. Thus, the single cloud model can still reproduce the absolute brightness 
and temporal variability seen in the low spectral resolution, visible EPOXI data, as 
shown in Table II and in the filters presented in Fig. 7 (the two cloud model can 
reproduce the visible data, and is not shown for clarity). These results emphasize the 
crucial role that spectra play in disentangling the effects of clouds; the EPOXI visible, 
broadband data can be fit with a single cloud category while four cloud varieties were 
required to fit the NIR data. It may also be possible to tweak the optical depths of 
the clouds in the two cloud or the four cloud model to improve the fit to the NIR 
spectra while remaining consistent with MODIS data. 


19 



5.2.2 Sensitivity to Atmospheric and Surface Resolution 

The atmospheric resolution determines the scale at which temperature and gas mixing 
ratio profiles are resolved while the surface resolution determines the scale at which 
surface features are resolved. To test our sensitivity to these parameters, we ran a 
model with a single atmospheric pixel as well as a model with 48 surface pixels (models 
“f” and “g” in Table II, respectively). Using a single atmospheric pixel amounts 
to assuming there is no spatial variability in the temperature and composition of 
Earth’s atmosphere. Earth’s surface exhibits large variations in both water vapor 
mixing ratios and surface temperatures, making a single atmospheric pixel a poor 
choice for a model that aims to simulate these variations. For example, a model with 
a single atmospheric pixel would incorrectly produce polar and equatorial thermal IR 
spectra that are nearly identical, even though, in reality, these regions can differ in 
temperature by over 100K. In general, though, the model with a single atmospheric 
pixel can still reproduce the visible EPOXI data reasonably well since the visible data 
is relatively insensitive to atmospheric/surface temperature as well as water vapor 
distribution. Regarding surface resolution, the coarse resolution of the 48 pixel model 
(surface pixels are several thousand kilometers in size) leads to a poor reproduction of 
the longitudinal variations in brightness in the lightcurves at all wavelengths, which is 
shown in Fig. 7. The visible disk of Earth in this low resolution model is dominated 
by r^j 10 pixels, which leads to unrealistic periodicities in the lightcurves as relatively 
bright surface pixels rotate into and out of view. 

5.3 Comparison to Reflectance Models 

To mimic previously published reflectance models (e.g., Williams and Gaidos, 2008), 
we removed atmospheric absorption and scattering from our model and replaced the 
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Mie scattering clouds in our model with a Lambertian surface with a visible reflectance 
of 0.60, which is a typical value assumed in reflectance models (model “h” in Table 
II). By definition, a reflectance model cannot reproduce spectral data, like the NIR 
EPOXI data, because such models explicitly ignore atmospheric absorption. The 24- 
hour average radiance is shown for the reflectance model, the EPOXI data, and our 
standard model in Fig. 9. The reflectance model clearly struggles to reproduce the 
radiance data (quantified in Table II), especially in the 350 nm and 950 nm Liters, 
which are strongly affected by extinction due to Rayleigh scattering and water vapor, 
respectively. However, the reflectance model can reproduce the relative variations in 
brightness in the EPOXI lightcurves, as shown in Table II and in Fig. 7. 

6 Discussion 

Earth is a complex system and, as a result, any model which aims to simulate Earth’s 
appearance to a distant observer should reflect this complexity. Our spectral Earth 
model aims to accurately simulate Earth’s disk-integrated spectrum at any arbitrary 
viewing geometry and wavelength, necessitating a comprehensive and rigorous treat- 
ment of a large number of physical processes (e.g., ocean glint, realistic cloud scat- 
tering, vertically and spatially resolved temperature and gas mixing ratio profiles). 
However, and as we have shown, this does not necessarily imply that more simpli- 
fied models cannot reproduce specific details of Earth’s appearance (e.g., ocean glint 
in Williams and Gaidos (2008)). In either case, any model of Earth’s appearance 
should be validated against observational data, ensuring its accuracy as a predictive 
or interpretive tool. 

Our new model is capable of reproducing the time variable color and absolute 
brightness of Earth, as observed in the visible and NIR EPOXI data, and can do 
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so for multiple dates of EP0X1 observations. Typical RMS errors for the model 
lightcurves are within 3-4%, and the 24 hour average visible radiance for the model 
matches the EPOXI observations to within calibration uncertainties. The model also 
simultaneously provides a good fit to the shape and absolute brightness of the high 
resolution AIRS mid-IR observations. 

Although our fit to the lightcurves is generally good (Fig. 3), the largest root- 
mean-square errors are seen for the 750 nm and the 850 nm filters, indicating that our 
poorest match to the shapes of the EPOXI lightcurves occurs at these wavelengths. 
These filters are relatively clear of atmospheric absorbers and are largely unaffected 
by the strong Rayleigh scattering seen in bluer filters. As a result, these filters are 
the most sensitive to the surface, and mismatches in these filters may indicate that 
more than five surface types are needed to better reproduce the EPOXI data at these 
wavelengths. 

Even though our model reproduces the 24 hour average visible radiance of Earth 
to within instrument uncertainties (Fig. 4), discrepancies at short wavelengths are 
typically in the sense of the data being brighter than the model, suggesting either a 
systematic calibration error or residual minor defects in the model. The difference 
in the 24 hour average brightness between the data and the model is largest in the 
350 nm and 450 nm filters. A small fraction of the light incident on an ocean surface 
actually enters the water and is scattered back out (Cox and Munk, 1954), which 
is an effect not accounted for in our model. As water is most transparent in the 
350 nm and 450 nm filters, including this ocean “volume scattering” behavior in our 
simulations could improve our fits to the Earth’s radiance. 

Notable disagreements between the model and data in the NIR spectra occurred 
near 1.4 /mi, 1.9 /mi, and 4.1 /mi. The 1.4 /mi and 1.9 /mi discrepancies occur 
at the base of water vapor absorption features, indicating problems with the input 
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MODIS water vapor distribution data, the HITRAN line lists, and/or the vertical 
placement of clonds in our model. The height of clouds in the atmosphere controls 
the column depth of water vapor that is available to absorb radiation incident on the 
top of the atmosphere. For this reason, low clouds allow a longer column through the 
atmosphere and more absorption by water vapor in the spectrum than high clouds. 
MODIS data does not contain a complete description of the full three-dimensional 
distribution of clouds in Earth’s atmosphere and we must estimate standard altitudes 
for our liquid and ice clouds. While the altitudes that we assumed seem to offer a 
reasonably good fit overall to the NIR data, these fits could be improved by allowing 
our clouds to have a varying vertical distribution that is determined by CloudSat 
data (Stephens et ah, 2002). 

The disagreement near 4.1 //rn is in the wing of a CO 2 absorption feature. Efforts 
to fit this shape by altering cloud coverage and thickness were unsuccessful. N 2 0 
has a weak absorption feature between 4.0-4. 1 fin 1 , but altering atmospheric N 2 0 
levels also failed to reproduce the observed shape. S0 2 absorbs in this region, but 
the strength of this feature is even less than the N 2 0 feature. It is possible that 
we are missing a trace gas which absorbs in this region, or that our C0 2 linelist is 
incomplete, although a test using the more recent HITRAN 2008 database (Rothman 
et ah, 2009) did not improve our fit in this region. 

The original version of the model (Tinetti et al., 2006a, b) failed to reproduce 
the observed lightcurves primarily because this version of the model required a less 
rigorous parametrization of clouds to reproduce limited Earth-observing data sets. 
Most importantly, this earlier version of the model assumed that all ice clouds were 
quite thin with an extinction optical depth of order unity, which is true for only ~ 10 % 
of all ice clouds in the MODIS data. The model presented in this work has a much 
improved treatment of spatially and temporally varying clouds, and parametrizes 
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them based solely on input data collected from the MODIS instruments. Our cloud 
parametrization technique is rigorous and versatile, allowing us to reproduce all three 
sets of EPOXI observations, which span almost three months in time, without needing 
to tune model parameters to each dataset. 

In summary, our model is designed to be comprehensive and versatile enough to 
model the Earth’s appearance over a very large wavelength range and at arbitrary 
viewing angle and phase. The level of model complexity required to simultaneously 
simulate Earth’s spectrum over a large wavelength range, where different physical 
processes dominate, may at first appear daunting. Especially if considering the re- 
verse problem of retrieval of the correct planetary characteristics from a limited data 
set. However the more optimistic view is that in cases where a more comprehensive 
model is required to accurately fit Earth data, this indicates that the data contained 
enough information to allow us to discriminate the more complex environmental char- 
acteristics from more simplistic models. This would be a desirable circumstance when 
attempting to learn about extrasolar planet environments from observations. Addi- 
tionally, in the process of fitting the EPOXI data we have been able to quantify when 
model complexity is and isn’t required for a particular application or wavelength 
range subset. 

For example, modeling moderate resolution NIR spectra does require multiple cat- 
egories of clouds, providing cloud altitude, phase, and optical thickness resolution. As 
demonstrated in Fig. 8, a single cloud category produces a disk integrated spectrum 
that is too bright in the NIR continuum and underestimates water vapor absorption 
in some regions (near 1.4 fin i) while overestimating water vapor absorption at other 
wavelength regions (near 1.1 //m) . The lack of ice clouds leads to discrepancies near 
1.5 /irn where ice particles absorb. A model with two cloud categories reproduces the 
spectral data more accurately than the single cloud model, but struggles with the 
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shape of the 1.1 /mi and 1.4 n m water features. Residuals for the two cloud model 
can be over 40% larger than the residuals for the four cloud model in these regions. 
The four cloud model therefore appears to be optimal for simultaneously fitting both 
the visible and NIR spectral regions, and would be most useful for studying the de- 
tectability of Earth’s globally averaged characteristics for TPF-like designs that span 
both the visible and NIR. 

The absolute brightness and temporal variability of the EPOXI lightcurves can 
be reproduced by models without a large number of cloud categories due to the 
broadband nature of these data and the fact that observations at these wavelengths are 
relatively insensitive to atmospheric absorption. In essence, the broadband lightcurves 
provide evidence for white, highly reflective structures that vary in time on the planet, 
and observations in the 950 nm filter demonstrate an absorption feature from water 
vapor. The higher spectral resolution NIR data provide information regarding the 
phase and vertical distribution of these structures. Furthermore, insofar as both 
liquid water and ice clouds are required to reproduce the observations near 1.5 /mi, 
the broadband data and moderate resolution spectra demonstrate that water is found 
in the atmosphere as vapor, liquid, and ice. 

Sensitivity tests indicate that high atmospheric resolution is not needed to repro- 
duce the EPOXI visible photometric data. This is not surprising as at visible the 
temperature structure and distribution of trace gases within the atmosphere should 
have only small effects on the lightcurves. Even data in the 950 nm filter, which con- 
tains a large water vapor absorption feature, can still be fit due to variations in the 
brightness of the continuum outside the absorption feature and the fact that clouds 
control the column depth of water vapor that is available to absorb radiation. Earth 
exhibits large variations in both water vapor mixing ratios and surface temperatures, 
indicating that models with low atmospheric resolution are poor choices for modeling 
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high resolution spectral data or mid-IR data, especially if the model aims to generate 
observations for arbitrary viewing geometries (e.g, polar versus equatorial views). 

Reflectance models that ignore scattering and absorption in the atmosphere and 
which treat clouds as Lambertian reflectors cannot reproduce the 24 hour averaged 
brightness of Earth. The short wavelength Liters have lower reflectivity than the 
data because Rayleigh scattering has been ignored while longer wavelength filters 
{e.g., the 950 nm filter, which contains a strong water feature) show enhanced re- 
flectivity because atmospheric absorption has been ignored. Reflectance models can, 
however, reproduce the shapes of the EPOXI lightcurves as these models are designed 
to reproduce relative brightness variations due to structures {e.g., clouds, continents) 
rotating into and out of view. While simple, computationally inexpensive models, like 
reflectance models, may be useful as retrieval tools in scenarios where observational 
data are limited or of poor quality, these models are not optimal when compared 
to more rigorous and comprehensive spectral models for applications which require 
accurate predictions. 

Our validated model has a variety of applications. In Robinson et al. (2010), the 
model was used to demonstrate that surface oceans on Earthlike extrasolar planets 
may be detectable, even after considering the confusing effects of forward scattering 
from clouds. The Tinetti et al. (2006a, b) model could not perform such a study as 
it did not simulate specular reflection from ocean surfaces, and reflectance models, 
while excellent at reproducing glint, could not address the important issue of phase 
dependent aerosol scattering when investigating the detection of ocean glint in the 
presence of clouds (Williams and Gaidos, 2008). In Cowan et al. (2011), the model 
presented here was used to simulate observations of a distant Earth to validate and 
better understand a retrieval method for exoplanets. 

Our comprehensive model is uniquely capable of investigating a variety of Earth’s 
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traits over wavelength ranges, synoptic views, and vantage points that are unavailable 
to Earth observing spacecraft or satellites. For example, future applications could 
include model generated disk integrated, ultraviolet, visible or IR spectra of Earth 
for a variety of sub-observer points and phases. These simulated datasets could be 
used, among other applications, to investigate the wavelength dependent effect of 
clouds on our ability to measure thermal radiation from the surface, or to simulate 
Earth as seen from a lunar vantage point over a complete lunar orbit. We could, also, 
investigate the effects of an unresolved, airless satellite on the spectrum of its host 
by pairing our simulations with a model of the spectrum of the Moon. Finally, our 
Earth model could also be used to generate synthetic observations for “blind” tests 
of retrieval models, where other individuals or teams attempt to retrieve planetary 
characteristics without knowing the input to our Earth model ( e.g ., season, viewing 
geometry, phase). In general, simulated data can be used to test techniques aimed at 
characterizing habitable planets that may be employed by TPF-class missions. 

7 Conclusion 

We have developed and validated a three-dimensional spectral model of Earth that 
is capable of reproducing the temporal variability and absolute brightness of obser- 
vations in both the visible, NIR, and MIR. Earth’s visible lightcurves are strongly 
dependent on cloud spatial distribution and reflectivity. To simultaneously reproduce 
the EPOXI visible and NIR observations of Earth we used four categories of water 
clouds defined using data from the MODIS instruments. Smaller amounts of cloud 
categories can reproduce the visible EPOXI data (which is at a very low spectral res- 
olution, R ~ 6), but cannot reproduce the moderate resolution (R ~ 500) NIR data. 
Our simulations are relatively insensitive to variations in the surface resolution and 
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can also reproduce the EPOXI lightcurves at a very low atmospheric resolution. The 
model can now be used as a forward model to explore the detectability of planetary 
characteristics by generating synthetic observations of Earth from the far ultraviolet 
to the far IR at a variety of spectral, spatial and temporal resolutions. 
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10 Tables 


Tabic I: Summary of Trace Gas Input Data 


Species 

Instrument 

Valid Range [Pa] 

Typical Mass Mixing Ratio 
Surface Tropopause Stratopause 

ch 4 

Aura/TES 

1 x 10 5 - 5 x 10 2 

1 x KT 6 

1 x KT 6 

2 x KT 7 

CO 

Aura/MLS 

1 x 10 4 - 1 x 10 2 

IQ- 8 - 10“ 7 

0 

1 

OO 

1 

0 

1 

-D 

1 x 10" 7 

h 2 o 

AIRS 

1 x 10 5 - 1 x 10 4 

O 

1 

CO 

1 

0 

1 

to 

3 x KT 6 

3 x KT 6 

h 2 o 

Aura/MLS 

3 x 10 4 - 2 x KT 1 




n 2 o 

Aura/MLS 

1 x 10 4 - 1 x 10 1 

5 x 1(T 7 

4 x KT 7 

2 x 1(T 8 

0 3 

Aura/MLS 

2 x 10 4 - 2 x 10° 

IQ" 8 - IQ" 7 

O 

1 

-4 

1 

O 

1 

05 

5 x KT 6 
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Table II: Sensitivity Test Results 


Model 

a 

b 

c 

d 

e 

f 

g 

h 

Observation Month 

March 

May 

June 


March 

March 

March 

March 

March 

Atm. Res. [pixels] 

48 

48 

48 



48 

48 


1 


48 

n 

/a 

Srf. Res. [pixels] 


192 

192 

192 


>3,000 

192 

192 


48 

192 

Cloud Cat 




1 


1 


1 



3 


1 


1 



1 

1 (Lambert) 

Filter Center 

nm 


24-hr Avg. Rad. Error 1 

[%] 

: Lightcurve RMS Error 2 

[%] 


350 



6.7 

1.4 

8.1 

3.4 

5.5 

3.3 

3.1 

: 4.5 

4.8 

1.7 

6.5 

1.4 

1.0 

1.6 

31 

3.1 

450 



7.0 

2.0 

8.1 

3.4 

7.1 

3.5 


1.3 

: 5.3 

3.2 

2.3 

6.7 

1.9 

0.9 

2.1 

13 

1.6 

550 



2.1 

2.3 

4.0 

3.8 

2.9 

3.8 

9.8 

: 5.1 

3.6 

2.6 

0.8 

2.2 

4.8 

2.5 

19 

1.4 

650 



4.2 

2.0 

5.8 

3.7 

5.3 

3.7 

7.6 

: 4.2 

2.8 

2.3 

3.0 

1.9 

2.5 

2.6 

27 

1.7 

750 



0.5 

4.3 

0.5 

5.7 

0.5 

5.8 

5.7 : 5.1 

8.7 

3.8 

0.5 

4.1 

6.6 

2.9 

34 

2.3 

850 



0.2 

5.5 

1.3 

7.2 

0.3 

7.2 

2.6 

: 6.6 

8.4 

4.6 

0.4 

5.4 

5.2 

3.7 

27 

2.7 

950 



5.7 

2.1 

6.3 

4.2 

3.3 

4.0 

13.5: 5.1 

7.2 

2.8 

7.3 

2.2 

0.1 

2.3 

56 

3.0 


1 percent error for 24 hour average radiance, a measure of ability to reproduce visible radiance of Earth; compare to ~ 10% instrument 
absolute calibration uncertainty 

2 root-mean-square error for normalized model lightcurves, a measure of ability to reproduce lightcurve shape 



11 Figures 

Figure Legends 

FIG. 1. A true color image of Earth taken from the EPOXI data set (left) and 
from our model (right). Using date specific cloud coverage and optical depth data 
allows us to match cloud features in the model to cloud features in the data. 

FIG. 2. Comparison between the Tinetti et al. (2006a, b) model (old model, gray) 
and our new model (black). A measure of Earth’s reflectivity, taken as 7 r times the 
disk integrated radiance (in W / m 2 / fim / sr) divided by the Solar flux at 1 A.U (in 
W/m 2 /iJ,m), is shown for the planet viewed at full phase, gibbous phase, quadrature 
(he., half illuminated), and crescent phase, or, alternatively, phase angles of 0°, 45°, 
90°, and 135°, respectively. Both models use realistic cloud cover, and the data for the 
old model is taken from the left column of Figure 6 in Tinetti et al. (2006b). EPOXI 
observations taken at phase angles of 57.7° and 76.6° (dashed gray) and our model 
of the observations (dashed black) are shown in the gibbous and quadrature cases, 
demonstrating that our model correctly reproduces the brightness and spectral shape 
of the data at these phases. In general, the Earth model from Tinetti et al. (2006a, b) 
is about 100% to 400% too bright, and is too blue at phases near quadrature and 
crescent. Earth views generated by the Earth and Moon Viewer, first implemented 
by J. Walker (http://www.fourmilab.ch/cgi-bin/Earth/). 

FIG. 3. Lightcurves of data (solid) and our new and improved model (dashed) 
for EPOXI observations from March (left), May (center), and June (right). March 
observations begin at 2008-Mar-18 18:18 UT, May observations begin at 2008-May- 
28 20:05 UT, and June observations begin at 2008-Jun-4 16:57 UT. The filter center 
wavelength is noted in the central column. Model values and data have been normal- 
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ized to their respective 24 hour averages. The shaded region in the central column 
marks a Lunar transit of Earth’s visible disk, which is an effect not included in our 
spectral model. The vertical gray line indicates where the observations begin and, 
24 hours later, end. The discontinuity here is a real effect due to time varying clond 
structure and is of order 2-3% in March. The discontinuity tends to be smaller in the 
May observations and larger in the June observations. The model generally repro- 
duces the scale and sense of these discontinuities. 

FIG. 4. Comparison of the 24 hour averaged signal for the model (dashed) with 
the EPOXI data (solid) for the March (upper) and June (lower) dates of observa- 
tion, demonstrating our fit to the data on an absolute scale. Note that the June 
observations are overall dimmer than the March observations due to Earth phase. 
The largest discrepancies are typically in the 450 nm filter and are ~ 8% for both 
observations, within the 10% absolute error in the HRI calibration (Klaasen et ah, 
2008). The average spectrum of the May observations is similar to that of the June 
observations (he., within a few percent) and were omitted for clarity. 

FIG. 5. Near-infrared spectral comparison of the model (dashed) with EPOXI 
data (solid) for a variety of observations, demonstrating the ability of the model 
to reproduce moderate spectral resolution observations on timescales from hours to 
months. Note the different scales used for the y-axes on the left and right sides of 
the spectral plots and also note that the model has not been scaled to match the 
data. Date indicators and sub-observer longitudes are given at the top of each plot. 
Prominent absorption features have been labeled in the upper-left plot. May (middle- 
right) and June (lower-right) observations are dimmer in reflected light due to Earth 
phase. Residuals for the data- model comparison are shown below each plot and are 
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typically less than about 15%. Stronger disagreements tend to occur between 1.1- 
1.2 jum (water), between 1.35-1.5 /mi (water), and between 3.9-4.15 /mi (N 2 0, wing 
of C0 2 feature). Extremely low signal levels and instrument artifacts lead to large 
residuals in the 2.5-3.25 /mi range. Instrument calibration uncertainties are typically 
10%, and tend to increase below 2.0 /mi and above 4.3 /mi (Klaasen et ah, 2008). 
Earth views generated by the Earth and Moon Viewer. 

FIG. 6. Mid-infrared, 24 hour average spectra of Earth from our March model 
(dashed) and as generated from AIRS observations (solid) (Hearty et ah, 2009). In 
general, the agreement is quite good, with residuals (lower panel) being typically 
~ 7%. Large gaps are regions where the instrument does not return data. Note 
that the AIRS instrument does not record full disk observations, so disk integrated 
observations are generated from multiple scenes recorded over a 24 hour period. Thus, 
this technique cannot achieve the same time resolution as the EPOXI data set. 

FIG. 7. Comparison between the EPOXI data and a variety of models considered 
in this work through a subset of the EPOXI Liters for the March set of observations. 
Filter center wavelength is noted on each plot. The details of the models shown 
are outlined in Table II. Model “a”: standard model; model “d”: 2006 version of 
model Tinetti et al. (2006a, b); model “e”: single cloud category; model “f” : single 
atmospheric pixel; model “g”: 48 surface pixels; model “h”: reflectance model. Filters 
were selected to demonstrate the effects of Rayleigh scattering (350 nm) and water 
absorption (950 nm). The 650 nm Liter is relatively free of atmospheric extinction. 

FIG. 8. Comparison between the EPOXI data (solid), our standard model 
(dashed), a model run with a single cloud category (dotted), and a model run with 
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two cloud categories (dot-dashed) for a view over the Pacific Ocean on 2008-Mar-18 
UT. Note the different scales used for the y-axes on the left and right sides of the 
spectral plots. While the single cloud model and the two cloud model can reproduce 
the visible EPOX1 lightcurves, they cannot reproduce the NIR data. 

FIG. 9. Comparison of the 24 hour averaged signal (top) for the EPOXI data 
(solid), our standard model (dashed), and a model where atmospheric absorption and 
scattering has been removed and clouds have been treated as a Lambertian surface 
with an albedo of 0.60 (dotted). The data and models have been converted to a 
measure of reflectance (bottom) in the same fashion as in Fig. 2. The effects of 
ignoring Rayleigh scattering can be seen in the shortest wavelength Liters while the 
lack of atmospheric absorption is especially apparent in the 950 nm filter, which 
includes a strong water absorption feature. Data and models are all for the March 
observations. 
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